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ABSTRACT 
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reported. Measurements were made downstream using a small 
calibrated pneumatic probe, with simultaneous measurements 
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and temperature. In addition, a program is presented to 
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I. INTRODUCTION 



The study reported here is part of an ongoing program 
at the Naval Postgraduate School to determine the perform- 
ance of a transonic compressor stage. The overall goal of 
the program is to develop the knowledge necessary to design 
compressors with multiple transonic stages which are effi- 
cient, light in weight, and free of flutter problems. It 
has long been appreciated that increasing the wheel speed 
and relative air velocities in axial compressors would re- 
sult in fewer stages for a given energy addition. However, 
in order to achieve high overall efficiencies in multiple 
axial stages, the ability to predict the behaviour of a 
single stage correctly must be shown. The purpose of this 
study was to measure the blade -element performance of a 
transonic compressor rotor operating at high subsonic speeds 
and to compare the measurements with the performance calcu- 
lated according to the method given in Reference 1, 

The transonic compressor and drive turbine were de- 
signed by Dr. M. K, Vavra. The program is supported by the 
Naval Air System Command, Code 31 0» through the office of 
Dr. K. J, Mueller, 

To determine the blade-element performance of a com- 
pressor blade row, it is necessary to evaluate radial velocr-- 
ity and angle distributions into and out of the blade row. 

In a transonic compressor, at or near design speed, the 
unsteady r^ture of the flows with large pressure amplitudes 



14 



complicates the problem. In addition, the possibility that 
wakes might induce vibrations in the highly stressed blades 
prohibit the use of probes upstream of the rotor. Measure- 
ments can be made downstream of the rotor; however, since 
the flow area is small, boundary effects on probe measure- 
ments can be significant. 

In order to carry out the present study, a small com- 
bination temperature -pneumatic probe was designed by Dodge 
(Ref, 2) specifically to be used immediately downstream of 
the rotor. In previous work, the flow upstream of the rotor 
was determined at different flow rates (Ref, 3 )* It was 
left to the present work to represent the inlet flow field 
analytically, to calibrate the combination probe for bound- 
ary effects, and both to measure and to calculate the rotor 
blade element performance. 

In this report the transonic compressor and test rig 
are first described in Section II. The method of approach 
and program of measurements are described in Section III. 

The results of the measurements are then described and 
discussed in Section IV, and conclusions are given in Sec- 
tion V. In order to maintain a clear presentation of the 
relevant results, the details of the calculations and of 
the calibration measurements are given in separate Appendices, 
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II. THE TRANSONIC COMPRESSOR 



A. DESCRIPTION 

A drawing of the compressor test rig is shown in Fig. 1. 
The test rig consists of an air-turbine drive unit and an 
induction section containing a filter, a throttle, a set- 
tling chamber and a flow measuring nozzle. The air turbine 
is designed to deliver 450 horsepower at 30»000 RPM. The 
relative blade tip Mach number of the rotor at design point 
is 1.5* The compressor flow rate is controlled by an elec- 
tric-hydraulic rotating throttle plate. Power is provided 
by a compressed air supply from an Allis-Chalmers multi- 
stage axial compressor. A complete description of the 
Turbopropulsion Laboratory and test facilities is contained 
in Ref. 4. A drawing of the transonic compressor is shown 
in Fig. 2. 

The blade design for the transonic compressor rotor was 
derived from circular-arc elements wrapped on six conical 
surfaces. The outer surface (at the case wall) is a cylin- 
der. The blading geometry is given in Fig. 3. The blade 
pressure side is flat and the suction side is a circular- 
arc. The solidity is fairly constant with a value of 1.35 
at the mid-streamline. The outer diameter of the rotor 
blade is 11 inches with a channel height of 1.92 inches 
downstream and a radius ratio at the inlet face of 0.5. 
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B. INSTRUMENTATION 

The transonic compressor is equipped with stationary 
probes to provide routine performance data. Triple redun- 
dancy is provided in the measurement of the power; either 
the mass flow rate and temperature drop through the drive 
turbine, or the mass flow rate and temperature increase 
through the compressor, or the speed and the measurement of 
torque can be used. The torque is measured by strain gauges 
on flexures against which the stator assembly is free to 
rotate. Flow straighteners before and after the stage en- 
sure that the stator and rotor torques are equal. The 
method based on torque and speed is the most accurate since 
it is independent of flow profiles. Static pressure taps 
were installed on both outer and inner walls of the com- 
pressor casing. The rotative speed of the test rotor was 
determined by a magnetic flux cutter and electronic counter. 
The fixed instrumentation will not be described further. 

It was used in the present tests only to provide the refer- 
ence conditions for the survey measurements. The compressor 
cylindrical casing was modified to accommodate a variety 
of probes and sensors which may be required in the program 
of steady state and unsteady flow pressure measurements. 

The casing can be rotated to position probes in the periph- 
eral direction. 

The rotor exit survey measurements reported here were 
obtained using a combination probe (Fig. 4) designed by 
Lt. F. J. Dodge (Ref. 2). The combination probe, following 
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calibration, was used to determine total pressure rise, to- 
tal temperature rise, and flow velocity in magnitude and 
direction. In addition, static pressure was measured using 
the wedge static probe shown in Fig. 5f which was based on 
a NASA design provided by Mr. Lloyd N. Krause of NASA Lewis 
Research Center. Each probe had two "static pressure" sen- 
sors which, if the probes were rotated until the indicated 
static pressures were equal, allowed the probes to be aligned 
with the flow direction. It was determined that better 
accuracy could be attained, however, by aligning both probes 
to the flow direction sensed by the combination probe. This 
was because the sensors on the wedge probe were separated 
in the radial direction. A chromel-constantan fine-wire 
thermocouple in the combination probe was used in conjunction 
with a similar sensor mounted upstream of the rotor to 
measure, differentially, the total temperature rise. 

Radial surveys were taken at a fixed peripheral case 
angle by manually positioning the probe in the flow using a 
United Sensor traverse unit, manually rotating the combina- 
tion probe to balance indicated static pressures, and ad- 
justing the wedge probe to the same angle. All pressures 
were measured by a calibrated transducer in a Scanivalve 
arrangement. (A capability exists to remotely control probe 
movement and balancing through the use of a programmable 
microprocessor and a probe traversing mechanism designed by 
Lt, D. D, Patton (Ref, 5). Remote operations will be needed 
as higher operational speeds are developed.) 
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C. DATA ACQUISITION AND REDUCTION 

Pressiires, temperatures, rotational speed and other data 
were scanned and punched onto paper tape hy a B & F Model 
24133 solid state data acquisition system. A complete de- 
scription of this system is contained in Ref, 4, The 
punched tape data was processed using Hewlett-Packard Model 
983 OA programmable calculator data system. The data was 
transferred to a Model 9867 B mass memory disk unit for stor- 
age prior to data reduction. The data facilities are shown 
in Fig. 6, 

A program, designated "\W001" was written to transfer 
raw data from tapes to the mass memory. A second program 
"HAWKl" was written to calculate temperature, corrected 
non-dimensional velocity, turning angle, total pressure, 
and static pressiure. Two other programs, used in sequence 
("REPROl" and" "RA5C") completed the data reduction process. 
Plotting routines developed by Dr. R. P. Shreeve were 
adapted to plot distributions of Mach number, pressiare rise, 
temperature rise, flow angle, and loss coefficients, using 
the X-Y plotter. 

In reducing the probe measurements to obtain velocities 
and loss data, the thermocouple calibration reported in 
Ref. 2 was used. However, new pneumatic calibrations were 
performed and corrections were derived as described in 
Appendix 3, 
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III. METHOD OF APPROACH 



A . METHOD 

The usual approach to the design of axial compressors is 
described in Ref, 1. The method involves combining an axi- 
symmetric through-flow analysis, which prescribes the stream 
surfaces, with a two-dimensional cascade performance calcu- 
lation which determines the required blading geometry on 
those surfaces. The blading is generally layed out on coni- 
cal approximations of the prescribed stream siorfaces. 

When a new machine is built, the opportunity is there 
to compare the acttial blade element performance established 
by measurement with that which can be calculated by applying 
the design approach. The comparison can be used to derive 
corrections for use in future designs, and it is particularly 
useful to obtain data for new blade profiles. 

The transonic compressor investigated here was designed 
using the above approach, and the resulting blading is shown 
in Fig. 3. It is to be noted that the rotor blading is not 
typical in design. In order to establish the blade element 
performance of a rotor experimentally, measurements should 
be made of the flow into and out of the rotor on the stream 
sxirfaces on which the blading profiles are prescribed. How- 
ever, it is desirable not to operate at high wheel speeds 
with strong wakes from upstream probes, and so in the pre- 
sent work a method of prescribing the inlet flow field from 
measurements of only total flow rate and wheel speed was 
derived. The method is given in the following section. 
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The flow downstream of the rotor was determined using two 
calibrated probes. The calibration and corrections re- 
quired for the probes are described in Appendix B, From 
the rotor exit flow measvirements, the blade element loss 
coefficient and the local deviation angle were calculated. 

Finally, as described in Appendix A. the method of 
calculating deviation angle and loss coefficient given in 
Chapter 7 of Reference 1 was programmed for the Hewlett- 
Packard Model 983 OA calculator so that these quantities 
could be computed for four blade element surfaces at any 
given compressor flow rate and operating speed, 

b ', inlet flow field representation 

The flow field ahead of the rotor was measured in early 
tests using a calibrated pneumatic probe (Ref. 3). The 
measurements were compared with the results of a compress- 
ible finite difference calculation reported in Reference 6, 
and blockage factors were derived to relate the inviscid 
calculations to the measxired flow field. 

The method adopted here was based only on Anderson's 
measured data. As described in Appendix C a polynomial re- 
lation was found to correlate the measured velocity distri- 
butions away from the hub and tip regions. The flow function 
in this region was therefore calculated using Eq. C(3) and 
Eq. C(7). 
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C. USE OF FLOW PROBES 

Pressure probes are designed to exploit the distribution 
of pressure which occurs over a body immersed in a moving 
fluid. These pressure variations depend mainly on the flow 
velocity so that with a carefully selected choice of body 
shape, having multiple sensors^ a probe can be calibrated to 
determine flow velocity. The relationship between pressure 
and velocity must be established over a range of Mach num- 
bers. In addition to flow velocity in magnitude and direc- 
tion, the distributions of stagnation temperature, and total 
and static pressure are usually required. The probe designed 
by Dodge and reported in Ref, 2 is a small lightweight probe 
of aerodynamic design which allows these quantities to be 
determined at any position in the flow. 

In order to provide an independent measurement to com- 
pare with the -combination probe, the static pressure distri- 
bution was measured with a second probe, and static pressures 
were also recorded at hub and tip wall taps. The wedge 
static probe shown in Fig. 5 has the favorable characteristic 
that when it is aligned with the flow in the yaw plane (by 
balancing the pressures measured on the wedge surfaces) the 
indicated static pressure is relatively insensitive to angle 
variation in the pitch plane. As reported in Ref, 7 this 
type of probe has been used successfully up to a Mach number 
of 1.6. Fxirther, for measurements in the transonic regime 
the probe has favorable characteristics. If the leading 
edge is sharp and the wedge is small. Ref, 7 states that the 
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calibration is also independent of Mach mamber. This was 
found to be the case for the probe used in the present work. 

D. PROGRAM OF MEASUREMENTS 

The combination probe and the wedge static probe were 
mounted at the exit of the rotor with a peripheral separa- 
tion of 120 degrees, Foiar surveys were conducted at 50 
percent of design speed, at throttle conditions varying 
from near surge to maximum flow rate. Case angle for the 
first survey was 3^5 degrees. Subsequent svirveys were con- 
ducted at a case angle of 357 degrees because it was deter- 
mined that peripheral variations induced by the stator 
blading were felt least at this setting, RPM was the same 
(15»230t30) for all four runs. 

The procedure in each survey was to move the probes in- 
ward simultaneously, balance the combination probe to sense 
the flow angle, and then set the wedge probe to the measured 
angle, A data scan was taken, and the procedure was repeated 
at 0.1 inch intervals. Because of the hooked geometry of 
the wedge probe, measurements with this probe of static 
pressure could not be made closer than ,6 inches to the 
outer wall boundary. Static pressures were derived for the 
complete passage through the use of the combination probe 
calibration. Wall static pressure measurements at the hub 
and tip were also recorded. Yaw angle and linear displace- 
ment were recorded manually. All other data were automati- 
cally recorded by the data recording system. Pressures were 
recorded to an accuracy of .1 inches of water. 
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V. RESULTS AND DISCUSSION 



A. TIME-AVERAGED FLOW FIELD 

The distributions of total pressure rise, total temper- 
ature rise, Mach number and flow angle measured using the 
combination probe downstream of the rotor are shown in 
Figs. 7 1 8, 9, and 10. 

The total pressure rise was nearly constant with radius, 
whereas the temperature rise was indicated to be larger 
toward the tip. The larger indicated temperature rise near 
the case wall might be caused by higher temperatures in the 
incoming case wall boundary layer. Both the Mach number 
and the flow angle increase toward the hub, so that the 
through-flow velocity tends to remain constant. The results 
for the flow angle are qualitatively consistent; however, a 
large decrease can be observed in the magnitude of the flow 
angle for the open throttle case (Fig. 8). No definitive 
explanation was found, but as will be shown later this 
measurement leads to inconsistencies in the correlation of 
blade element performance. 

The velocity and pitch angle distributions, which were 
derived from the probe calibration and corrected for bound- 
ary effects (Appendix B), are shown in Fig. 11 and Fig. 12. 
The pitch angle shovis a smooth variation between the values 
of the v/all slopes at the hub and tip boundaries. It is 
shown in Appendix B that the pitch angle is determined to 
only l|-° accuracy, so that the variations in pitch angle 
indicated in Fig. 12 might be exaggerated. 
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A sensitive test of the accuracy of the combination probe 
measurements- was to examine the static pressure distribution 
inferred from the measured velocity and probe impact pres- 
sure measurement. Fig. 13 shows the distribution of static 
pressure calculated in this way. Also included in Fig. 13 
are the values of static pressure measured at the hub and 
tip wall static taps. It can be seen that agreement near 
the case wall was excellent. The comparison near the hub 
wall was not as good) however, very similar results were 
obtained from the wedge probe as shown in Fig. 14. While 
no readings were possible near the outer wall using the 
wedge probe, the two probes gave very similar results near 
the hub wall, and gave almost identical results for the sur- 
vey at open throttle. 

As shown in Appendix 3, the effect of flow boundaries on 
the readings of the combination probe is small. The magni- 
tude of the correction which was included in the above re- 
sults can be seen in Fig. 15. The correction is less than 
in velocity over the useful range of the survey. 

Velocity diagrams near the mean radius (0.55 of the 
channel height) are shown in the four parts of Fig, 16. 

B. DEVIATION ANGLES 

The design procedure for calculating the deviation 
angles for a given blade profile, with corrections for 
three-dimensional effects, is described in Appendix A, 

31ade profiles are known on four conical surfaces in the 
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flow field, on which the blading dimensions are those given 
in Fig, 3. Using relative flow angle from the velocity- 
diagram data (Fig. 16 for example), the flow deviation angle 
(S) can be calculated from the stagger angle {Y) and camber 
angle using the relation 




The results for the four probe surveys are shown in 
Fig, 17'in comparison with the results of calculations 
following Ref, 1. 

Departures from values predicted using Ref. 1 are quite 
large. However, the design procedure strictly applies only 
to the case where the operation is at the point of minimum 
loss for the given blade element, A closer correlation 
might be expected as the higher operating levels are achieved 
and closer to design conditions are attained. (Ref. 1 states 
that at transonic inlet Mach number levels, the point of 
minimum loss should constitute the desired design setting. ) 

It is noted that the through-flow velocity is considerably 
increased through the rotor at the design speed setting, 

and this would effect the deviation angle. No explanation 
was found for the indication of negative values of deviation 
angle near the hub, although a careful check was made on the 
accuracy of the probe yaw angle measurements and on the ac- 
curacy of the data reduction program. It must tentatively 
be concluded that there are strong three-dimensional effects 
present v/hich can lead to the indication (locally) of nega- 
tive deviation angles. More data must be obtained to investi- 
gate this speculation. 
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C . LOSSES 

Loss coefficients for the flow over the blade elements, 
assiming that the stream surfaces were the conical surfaces 
given in Fig, 3» were calculated from the probe measurements 
using the methods described in Appendix D. The calculated 
loss coefficient ) was consistent with the definition 
of loss coefficient used in Ref. 1. Of the two methods 
described in Appendix D, Method 1 uses both the stagnation 
pressure rise and stagnation temperature rise measurements 
of the probe. Method 2 does not use the stagnation temper- 
ature rise measurement, but uses a temperature rise calcu- 
lated from the velocity and angle measurements, through 
Euler's turbine equation. 

The results are shown in Fig, 18. Good agreement was 
obtained bet%veen the two methods of calculating the loss 
coefficient, however negative values were calculated for the 
inner stream surface. It should be noted that all measured 
data and all blade properties varied smoothly with radius, 
so that smooth curves between the individual data values 
could be drawn. The distributions are in qualitative agree- 
ment with earlier measurements, and it is noted that nega- 
tive values of losses were previously reported by Dodge 
(Ref, 2) from measiurements made with the same combination 
probe. 

A negative value of the loss coefficient can result only 
from a measurement of stagnation pressure rise which is too 
high, or a measurement of stagnation temperature rise which 
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is too low. This can be seen if Eq, D(13) is rewritten in 
the form 



_/ 

(V = 



[l * 



[l ^ \z] 



^2 

2 iV(Y-i; 



r 2, y/(y-i) 

- t - ^iJ 



< 1 - 



J2- 



¥ 



2V 



where*^is the ratio of total temperatures, ) ♦ 









Clearly, since Xu 2 a-nd Xu]_ are positive , is negative only 

Since the prediction of Method 1 and 
Method 2 agree, and because a favorable comparison with ear- 
lier unreported measurements of temperature rise was made, 
it is tentatively concluded that the impact pressure rise 
measurement is too large. A positive error in impact pres- 
sure could result from the pneumatic integration of the peri- 
odic unsteadiness of the flow from the rotor. All parts of 
the probe are the open ends of continuous tubing. It is 
known that pneumatic probes of different designs integrate 
differently in periodic flow conditions. An investigation 
of the pneumatic averaging characteristics of the probe is 
required in order to correct the measurements of losses. 

It should be noted that the velocity and pitch angle calcu- 
lated from the time-averaged probe readings should not be 
as sensitive to the averaging behaviour of the probe. This 
is because all parts are made from similar tubing and will 
average in a similar way. The loss coefficient, however, 
is very sensitive to errors in impact pressure measurement 
alone. 
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While the magnitude of the loss coefficients is clearly 
in error, it is of interest to attempt to correlate the 
measurements on the basis of diffusion factor as defined in 
Ref, 1, The result of this attempt is shown in Fig. I 9 . The 
diffusion factor (D) was calculated from the definition 

W2 R-iWui + Rp^'^u? 

D = 1 _ + ^ 2__2 

Cr (R^ + R2) W^ 

where W denotes the relative velocity, R is the radius, is 
the. solidity and subscripts 1 and 2 denote stations upstream 
and downstream of the rotor, is the peripheral component 

of the relative velocity. The "total loss parameter", 

( cD cos^g/^^^) is consistent with that shown in Fig. 203 
of Ref. 1 and was calculated from probe measurements. 

It can be seen in Fig. 19 that a reasonable correlation 
of the loss coefficient data on the basis of diffusion factor 
was obtained for three of the four sets of data. The data 
for full open throttle departed significantly from the data 
from three other test conditions. It should now be recalled 
that the flov/ angles measured for this test were also incon- 
sistent, More data is again needed, and a correction must 
be devised for unsteady flow effects on the probe measure- 
ments before a complete correlation can be obtained. 
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VI. CONCLUSIONS 



The attempt was made here to determine the blade element 
performance of an "impulse" rotor designed for transonic flow 
(design point) but operating at high subsonic (off -design) 
conditions. The results were preliminary, but served to 
define the problems involved in obtaining necessary design 
data from empirical correlations of measurements made in 
actual machines. 

The results obtained depended on the performance of a 
particular combination probe when used downstream of the 
rotor. The follov/ing conclusions were drawn concerning the 
probe measurements : 

(i) The effect of the presence of flow boundaries on the 
probe measurements was small, but corrections were derived 
and applied in the present work. Static pressure derived 
from the probe measurements, static pressure from wedge probe 
measurements and wall static pressures were in reasonable 
agreement, 

(ii) The temperature rise across the rotor measured by 
the probe, using the probe calibration, was consistent with 
the rise deduced from Euler's turbine equation. All indica- 
tions were that the probe provided an accurate measurement 
of the temperature rise. 

(iii) The impact pressure measured by the probe in the 
unsteady flow downstream of the rotor was probably too high. 
An investigation of the unsteady response of the tubing used 
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in the probe design is required in order to correct the 
data obtained so far. 

(iv) The flow velocity derived from the probe measure- 
ments was probably not as sensitive to unsteady flow effects 
since all the pressure sensors on the probe were constructed 
from the same tubing. The derived velocity depends mainly 
on differences between the pressures registered by the sen- 
sors. In steady flow the velocity was determined to better 
than 1.0^ accuracy, to nearly sonic speeds. 

(v) The pitch angle was determined by the same probe to 
within l§-°. 

(vi) The wedge probe provided a good measurement of 
static pressure in the compressor when rotated to a pre- 
measured flow angle. 

As a result of (iii), the loss coefficients determined 
from probe measvirements were too small, becoming negative 
over the blading near to the hub. A correction to the im- 
pact pressure measurements and more data are needed before 
a useful correlation of the losses can be obtained. However, 
it was found that the present results from three tests and 
four different blade-element surfaces correlated reasonably 
well in terms of the NASA diffusion factor. 

The deviation angles measured in the present tests de- 
parted significantly from the values predicted using NASA 
SP-36. No explanation was found for the negative values 
measured near to the hub, although the possibility exists 
that they could be caused by three-dimensional effects in the 
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flow through the rotor. More data are needed to confirm the 
trends in the deviation angles. Agreement with the predicted 
values was not expected since the blade profiles were unusual, 
and were operating well away from the design condition. 
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FIG. 1 THE COMPRESSOR TEST RIG 
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FIG. 2 TRANSONIC COMPRESSOR ASSEMBLY 
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FIG. 4 COMBINATION PROBE (SHOWN WITH THERMO- 
COUPLE REMOVED) 
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FIG. 5 WEDGE PROBE 
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FIG. 6 SCHEMATIC OF THE HEWLETT-PACKARD MODEL 983 OA CALCULATOR DATA REDUCTION SYSTEM 




I 



I 



t 



I 





40 



(IMn» JBNH uru 




41 



FUM msu (fiLPHR) 



« 









* I t 









( 

I 





42 
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FIG. 10 RESULTi OF nU»E SURVEYS 01 
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THE COMBINATION PROBE 
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NORMALIZED DISPLACEMENT (HUB TO TIP) 

FIG. 12 PITCH ANGLE DISTRIBUTION DOWNSTREAM OF THE ROTOR 
USING THE COMBINATION PROBE 
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NORMALIZED DI5PL0CEHENT (HUB TO TIP) 

FIG. 13 STATIC PRESSURE DISTRIBUTIONS DOWNSTREAM OF THE ROTOR 
DERIVED USING THE COMBINATION PROBE CALIBRATION AND 
IMPACT PRESSURE MEASUREMENT 
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NDRHHLI2D DI5PUCEHENT (HUB TD TIP) 

FIG. 14 STATIC PRESSURE DISTRIBUTIONS MEASURED DOWNSTREAM OF 
THE ROTOR USING THE WEDGE STATIC PROBE 
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t PRDfiE DIBHETERS FRDH OUTER NOLL BOUNORRY OE TRONSDRIC CDHPRE55DR 

PIG. 15 BOUNDARY EFFECT CORRECTION APPLIED TO THE COMBINATION PROBE RESULTS 
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{^nTDR entrance: plane 
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FIG. l6d VELOCITY TRIANGLES AT MIDSTREAM FOR RUN 52 POINT 46 
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FIG. 17 DEVIATION ANGLES FROM PROBE MEASUREMENTS AND CALCULATIONS 
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D- I RF-Ue I OKI RRCTUR 

FIG. 19 DIFFUSION FACDOR CORRELATION OF LOSS COEFFICIENTS 



APPENDIX A t BLADE ELEMENT PERFORMANCE PREDICTION 



Al. INTRODUCTION 

Flow through the blading of an axial-flow compressor is 
an extremely complicated three-dimensional process. As a 
result, simplified design models have been devised that per- 
mit the use of two-dimensional techniques. The approach is 
described fully in Reference 1. Involved is first the cal- 
culation of the "through-flow" assuming axial symmetry, then 
the calculation of the "blade-element flow" assuming approx- 
imately two-dimensional behaviour. The latter calculation, 
based on Chapters 6 and 7 of Reference 1, is the subject of 
this section. 

Five quantities shown in Figure Al determine the aero- 
dynamic behaviour of a cascade in ideal flow. Three of these 
define the geometry of the blade; they are, the shape of the 
blade, the orientation of the blades with respect to the cas- 
cade axis, and the solidity. The fourth quantity identifies 
the direction of flow ahead of the cascade and finally, the 
Mach number must be knov/n. When these quantities are speci- 
fied, the inviscid flow may be determined. 

Real cascade flows however, experience many other in- 
fluences including friction, t.urbulence, three-dimensional 
and unsteady effects, and boundary layer interaction. The 
procedure outlined in Reference 1 correlates cascade and com- 
pressor blade element performance data. Various correction 
factors are introduced which account empirically for the 
effects which can not be predicted analytically. The result 
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is a method of calculating the perfomnance of compressor 
blading. The method can be applied to both rotors and stators 
over a fairly wide range of Mach numbers. The method is 
limited however to the calculation of performance close to 
design point. 

NASA SP -36 has become a reference text in the axial 
compressor industry. A brief summary of the method is pre- 
sented here, followed by a computerized solution which uses 
the extensive correlations of data presented in this manual. 

The main thrust was to develop a program to calculate 
the blade element performance of a blade design of arbitrary 
geometry. The program would provide a very powerful design 
tool, and would provide the basis for a program to predict 
the performance of a transonic compressor rotor. Reference 6 
reports an extension of this was used with considerable 
success up to a Mach number of about 1.5* 

A2. METHOD OF BLADE ELEMENT PREDICTION 

A2,l Two-Dimensional Incidence Angle . Any cascade of 
aerodynamic blading behaves qualitatively as shown in Fig- 
ure A2 where A i is the "range" of performance and the mid- 
point of the curve is defined as the reference minimum loss 
incidence angle (i^g^ ). For blading with a given camber 
( 0 * ) , thickness, distribution, and maximum thickness ratio 
(t/c), the solidity (O’) and blade angle (V) complete the 
definition of the cascade geometry.' Given^ ^ and the inlet 
velocity the method in Reference 1 allows the determination 
minimum loss incidence angle, the deviation angle, and the 
loss coefficient. 



57 



The reference incidence angle for any blade element can 
be calculated from correlations of two-dimensional cascade 
measurements. Two-dimensional reference incidence angle for 
10^ thick blades is given as 

^2-D = A(l) 

where is a function of blade thickness distribution and 
(iQ)j^Q is the zero-camber incidence angle (function of ^ ^ 
and ^ ) , n is a slope factor given by 

n = ( i - iQ) / 2 _^ A(2) 

which is also a function of^^ and , 

In Ref. 1, (iq^IO given in Fig. 137 and n is given 
in Fig. 138. has the value 0.7 for double circular-arc 

aerofoil shapes. 

A2.2 Two-Dimensional Deviation Angle . As in the analy- 
sis of reference incidence angle, the correlation of blade- 
element deviation angle is expressed in terms of blade 
geometry by 

S 2-D “ A(3) 

where K^is a function of blade type (0.7 for double circular 
arc). zero-camber deviation angle (function 

of ^ ^ and(^), m is a function of^^, for the different basic 
camber distribution, and b is an exponent that is also a 
function of ^ 

Since deviation angle varies somewhat with incidence 
angle, a corrected deviation angle is defined as 

^ 2-D° " \ ^S^0°^10 ^^c ■ ^2-D^^ di ^2-D 
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where (d^/d^)2_Q is the slope of the two-dimensional varia- 
tion of deviation angle with incidence angle at reference 
incidence, and i is the actual incidence angle. 

Camber and turning angle are related by the equationj 

f" ?1 - ?2 5 + i 



or 

+ i - S a(5) 

Once the two-dimensional deviation and incidence angles have 
been established, the turning angle can be computed. 

A2.3 Corrections for Three-Dimensional Effects . The next 
step in the procedure is to correct reference deviation and 
reference incidence angles for three-dimensional effects. 

This is accomplished by defining the three-dimensional refer- 
ence angles 

,)deg A(6) 



^^ref^3-D " ^2-D ^ ^ " ^2-D' 



and 



^Sref^3-D " ^2-0 ^ ~ 

where (i - ^ 2 -^^ " ^2-D^ functions of Mach number. 

In summary, the desired turning angle is given by Eq, 

288 of Ref. 1. 






^^c ■ Sz-D^ ■ ^^c ■ ^2-D^ 



[l.O - ^z-d] " 



A(8) 



A2.4 Loss Coefficient . The loss coefficient for a blade 
element can be calculated from the correlation given in 
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Fig, 203 of Ref, 1, Here the "total loss parameter" is cor- 
related in terms of "Diffusion Factor", The calculation of 
losses is given in Appendix D, 

A3, POLYNOMIAL REPRESENTATION OF CORRELATIONS 

All values on the right side of Eq, A(8) may be deter- 
mined using data given in Figures in Ref. 1 once the geometry 
of the blade has been given and the inlet velocity calculated. 

Each of the unknown values of Eq, A(8) are functions of 
blading geometry and air inlet angle or blading geometry and 
Mach number, and are described by well-behaved families of 
curves. Fig, I68 of Ref. 1 has been chosen to illustrate 
the process of computerizing the solution, which involved re- 
ducing the data correlations given in figures to analytical 
expressions , 

Eleven values of air inlet angle ( ) were chosen for 

each solidity given in Fig, I68, The values were plotted 
as slope factor (m) vs. air inlet angle (^i) for constant 
solidity. The data points for each solidity were fit by a 
third order polynomial curve using the Hewlett-Packard 983OA 
calculator and X-Y plotter. For simplicity, a fifth order 
plotting routine was always used with the higher order co- 
efficients equal to zero if a lower order resulted in a more 
appropriate fit. The plot and polynomial curve fits for 
each solidity are shown in Fig, A3. The polynomial curve 
fit for each set of data points gave six soefficients for 
each curve. Thus for each m vs, plot, the six coeffi- 
cients Aq, A^, A^ were calculated for the polynomial 
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expression for m as a function of ^ ^ « 

m = Aq + + A^^^^ + A^^^^ + A^0^^ A(ll) 

at each value of solidity with 0^ converted to radians to 
minimize the error in higher order terms. The coefficients 
Aq, A^, ,A^ derived from the curve fits of Fig, A3 were 



plotted against each 


of twelve values of solidity as follows i 


Aq (I) vs. 


<r{i) 


, ,12points 


A^ (I) vs. 


<r(I) 


, . 12points 


• 

A^ (I) vs. 


(S' U) 


A(12) 

. ,12points 


This produced curves of the 


form 


■^0 ^0 ^ 


^ ‘ 




. 1 ^ 


1 12 


^ 13 ^ 1J\ ^ 15 


A^ - Sq 


aid- ^2 ^ 


+ 3.^ + 3.1^ 6 + a^ 


• 

• 

A = a 5 + 

Aj ap 


a^^d' + ag^d^ 


A(13) 

5 3 5 4 5 5 

+ + 3.i^^(^ + a^y 



which resulted in polynomial c\nr/e fits shown in Fig. A4, 

In this case the higher order coefficients A|^ and A^ are set 
equal to zero. When all the coefficients in Eq, A(13) are 
known, the slope factor m, for any given cT and can be 
calculated from: 

m = (a^ + a^<r + a.^(r + a^<r^ + a^q- + A(14) 

^ (a^^ + a^V + ^2^^^ 

:;+(aQ^ + a^^<j + a2^(5^ + ^3^(5^ 

This procedure resulted in a single equation (Eq. A(l4))for 
m in which the coefficients were analytic functions of 
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The coefficients are given in Tables A1 and A2. Due to the 
well-behaved nature of the curves, the slope factor m could 
be calculated to well within the accuracy of the original 
data contained in Ref, 1. 

All data necessary to predict performance for circular- 
arc blading was reduced to polynomials using a similar 
procedure. The polynomials were included in a computer 
program to calculate blade element performance for specified 
input geometry and flow conditions, 

A^, COMPUTER PROGRAM 

A4,1 Description , Program "NASA 36 " was written for the 
Hewlett-Packard 983 OA calculator to calculate the incidence 
and deviation angles for a rotor with meridional inlet ve- 
locity, The program follows the method of Ref, 1 through 
the equations given in Section A 2 , Where a value from a 
Figure in Ref, 1 is required, it is calculated using a poly- 
nomial approximation to the data as described in Section A3. 

The program is self-explanatory through the listing 
given in Table A3. The symbols are defined in Table A4, 

The program performs the blade element calculations for a 
specified compressor operating condition. The inputs are 
referred flow rate and speed, and temperatures upstream and 
downstream of the rotor. Through the method given in Appen- 
dix C, these inputs define the inlet velocity as a function 
of the radial position. It should be noted therefore that 
the present version of the program includes a description 
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of the transonic compressor inlet flow-field, while it 
allows changes in blading geometry to be input. 

A 4.2 Operation . Detailed information on the operating 
procedures for the HP-9830 calculator, mass memory and as- 
sociated equipment are found in Ref. 9 and Ref. 10 . The 
following is a step-by-step operating procedure for the 
program called "NASA36". 

1 , After loading "NASA36" into the Hewlett-Packard 
983OA calculator press "Run" and "Execute". 

2 , The calculator will then display a series of inputs 
required to determine the geometry of the blading. These 
are shown in Fig. 3 and arei 

a) Stagger angle 

b) Chord with sharp leading and trailing edges 

c ) Tip radius 

d) 'Hub radius 

e ) Diameter of cone through stacking radius 

f) Maximum Profile Thickness 

g) Radius of Profile Contour 

h) Blade camber angle 

i) Profile chord 

j) Angle of cone thru stacking radius 

3, The operating parameters of the compressor rotor 
are now required. These arei 

a) Revolutions per minute 

b) Total temperature upstream of rotor 

c ) Temperature downstream of rotor 
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d) Referred flow rate 

Outputs include deviation and incidence angles, Mach number 
before and after the rotor, and referred flow function (^), 
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Coefficients for Eq, A(ll) 
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SOLIDITY 

FIG. a4 variation OF OCEFFICIENTS OF POLYNOMIAL CURVES IN FIG. A3 



APPENDIX B: CALIBRATION AND CORRECTION OF PNEUMATIC PROBES 



Bl. FREE-JET APPARATUS 

The free- jet apparatus used to calibrate probes is 
shown in Fig. Bl, The device consisted of a 8 inch diameter 
pipe with the air supplied by the laboratory's central 
Allis-Chalraers Compressor. A ^.2 inch converging nozzle 
was used to increase velocity and attain Mach numbers up 
to about .9. To determine if the exiting velocity was uni- 
form and undisturbed by the probe mounting arrangement, 
surveys were conducted with a Prandtl probe over a range 
of 5 Mach numbers at , two axial locations. Surveys were 
also conducted with the same probe mounted and traversed 
from opposite sides of the jet. The results shown in 
Fig, B2 indicated that the velocity profiles were uniform 
over the core of the jet extending to a diameter of 3*8 
inches. As suggested in Ref. 11, the non-dimensional velo- 
city (X) appearing in Fig. 32 is defined by 



X = vA^ bU) 

, - ^ ^ o 

where ~ when T^ is in R. Following 

Ref. 11 further, the relationships for stagnation pressure, 

total temperature, and Mach number are, in terms of X given 




^ = (1 - X2)V(V-1) BO) 

^t 



( 
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V 
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M = 



X 



B(4) 



VTrt“ V^- 1 Vi - 



The values of X in Fig, B2 were calculated using the meas- 
urement of stagnation pressure and the assumption of atmos- 
pheric static pressure in Sq. B(3), 

B2. CALIBRATION OF THE COMBINATION PROBE 

Calibration of the pneumatic probe reported in Ref, 2 
can be represented by defining two functions of the non- 
dimensional velocity, X, and the pitch angle, (f) , If 
and ^2/3 indicated stagnation and indicated static pres- 

sure of the probe respectively, then define 



to be calibration functions when the probe is rotated to 



ured by the probe. 

The complete calibration procedure is given in detail 
in Ref, 3. The probe was recalibrated to confirm the re- 
sults given in Ref, 2 before investigating corrections for 
boundary effects, A range of Mach numbers (X = ,13 to ,4) 
and pitch angles (-20° to 20°) were used. With true static 
pressure taken to be atmospheric and true impact pressure 
that measured by a reference Prandtl probe, the value of 
X was calculated from Eq, B(3). As in Ref, 2, polynomial 




B(5) 



F^ = Fjzj (X,0) 




B(6) 



balance P 2 = P^. Pj^ is the fourth, off axis pressure meas- 
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represen'ta'bions were found for "the calibraijion equa"bions 
written in the form 

X = X(Fg,iZS) B(7) 

and 

= F^ (X,^) B(8) 

where , and are the ratios of probe pressure measure- 
ments defined in Eq. B(5) and Eq. B(6). It is noted that 
Eq, B(6) departs from the definition used by Dodge in Ref. 
2. The present method was found to give more accurate re- 
sults. The coefficients in Eq, B(7) and Eq, B(8) derived 
from the calibration test results are given in Table Bl. 

As a test of the analytical representation of the 
calibration, the calibration data v/ere used to calculate 
values of X and 0 to compare v/ith the actual values. Table 
32 shows the results of this comparison for the velocity. 
The error in X is generally less than Vfo and only reaches 
higher values at high positive or negative pitch angles 
( 16° or -8°) or at higher Mach numbers (.88), At this 
point the error increases to about 2^, 

Fig, B3 shows the error in the pitch angle in the 
range of pitch angle from -4° to +12°, The error is seen 
to be ll|-° and is not a simple function of either pitch 
angle or Mach number (or X). There appears to be a limit 
to the capability of resolving 0. The error was larger at 
high negative and positive values of pitch angles but these 
v/ere outside the intended range of use. 
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B3. BOUNDARY EFFECT ON COMBINATION PROBE MEASUREMENTS 



In application, the probe is not in plane flow away 
from all boundaries. The indicated velocity and pitch angle 
will depend upon several factors including boundary layer 
and stem immersion effects. The method of approach assumed 
here was that the probe always obeys the calibrations es- 
tablished and that errors are induced due to the perturbation 
of the probe itself on the flow near the boundary. Thus, 
in any flow the probe measures a velocity, X, and a pitch 
angle, which are the sum of the actual and an "induced” 
component. 

Three surveys were conducted in a 10 inch pipe up to 
the maximum attainable Mach number of ,55» The results of 
calculating velocity applying the free jet calibration given 
above are shown in Fig, 34. 

Introducing subscripts a(actual), i(induced), and 
m<measured), the velocity, X, may be represented by 



X = 
ra 



X + X, 




where the fractional 'error' 
the presence of the wall, is 




in X, 
given by 



Xa (1 B(9) 

which is induced by 



£ T = 



- K 

m a 



B(10) 



X is determined from the probe calibration in the usual 
m 

way. Xg^ is the required, correct velocity, and£"^ is the 
error to be established by calibration tests. 
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B(ll) 



The measxired velocity is corrected using 

X 

X = r— 21— 
a 1 + ^ 

t X 

when the error £ is known. 

In the same manner actual pitch angle may be repre- 
sented by 

0a " 0 " ^ 

In the 10 inch pipe tests, the actual velocity was 
calculated from the probe impact pressure and pipe wall 
static pressure using Eq. B(3) and the actual pitch angle 
was zero. Thus the "induced" error in X was calculated 
using Eq. 3(10) and the "induced" error in ^ was calculated 
using Eq. 3(12). The two errors, and , were then ex- 
amined as a function of distance from the pipe walls ex- 
pressed as a number of probe diameter?; (the characteristic 
dimension of the probe was taken to be 0.1 inches). It 
was found that the errors became negligible at about 9.5 
probe diameters from the boundary. In order to derive an 
expression of a suitable form, the errors were plotted as 
a function of (9.5-D)» where D is the displacement as a 
number of probe diameters. The velocity error was then 
approximated by a polynomial of the form 

^x = Ao + (9.5 - D) + A2 (9.5 - D)^ 3 ( 13 ) 

where D is measured from the wall through which the probe 
v/as inserted. The result is shown in Fig. 35. The error 
was less than 2 % at distances greater the 0.25 inches from 
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the wall. In order to extrapolate the correction to larger 
values of velocity, the error was represented as a function 
of X. The coefficients in Eq. B(13) were plotted as a func- 
tion of X, and the result is shown in Fig, B 6 . Aq, and • 
A 2 were seen to be almost linear functions of X and were 
therefore approximated by 



A 

A 

A 



0 

1 

2 



0 

II 


" V.95 


" V.95 


= 


+ c.x 


+ C^X 


0 


1 .95 


2^,95 
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+ D, X 


+ D^X 
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1 .95 


2\95 



2 

2 

2 



where X was the value of X at D = 9 . 5 . Values of 
X(corrected) minus X( uncorrected) verses immersion depth are 
shown in Fig. 37. Because of the aerodynamic design of the 
probe, the corrections necessary for boundary effects are 
small. 

Results for the induced pitch angle error are shown in 
Fig, 38. Since these errors were within the resolution 
capability of the probe ( 1 ^°), no correction was made. 



34, COMBINATION PROBE MEASUREMENTS IN STEADY SWIRLING FLOW 
Additional errors may be present when a pneumatic probe 
is immersed in a swirling flow,' Fig. B9 shows the apparatus 
which was constructed to investigate this possibility. The 
apparatus consisted of an annulus with a channel height of 
inches (10 inches in outer diameter). A total of 40 
static wall readings were available and a provision was 
made for probe surveys at different axial stations. The 



87 




1 



i 



swirl velocity was generated by a single stage of stator 
blading at the entrance to the annulus. Initial measure- 
ments suggested the possibility of boundary layer separation 
and reverse flow near the hub. In an effort to correct this 
problem, a plate was constructed to decrease the channel 
height at the exit, isolating the pressure at the hub wall 
from the surrounding atmosphere, and blocking backflow. 

Results are given in Fig, BIO which shows the static 
pressure derived from the probe calibration. Agreement 
is very close near the hub with a moderate departure near 
the tip at the highest Mach number, 

35. CALIBRATION OF THE WEDGE PROBE 

A type of probe which is reported to give a good indi- 
cation of static pressure over a wide range of Mach number 
is the v/edge static probe illustrated in Fig, 5 . The probe 
has a hatchet-shaped head with small pressure taps in each 
face of the wedge. In principle, if the leading edge is 
sharp and the wedge angle is small, the flow past the 
sensing orifice will be almost unperturbed, and the meas- 
ured pressure will be close to the static pressure in the 
undisttirbed flow. 

The static pressure may thus be determined by 

Ps = ( V3 - '’a > * Pa 

where it is assumed the probe is aligned in the flow with 
?2 equal, and is the reference atmospheric 

pressure , 
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The wedge probe shown in Fig. 5 was traversed in both 
the free jet and the swirling annular flows through the 
same range of Mach numbers reported for the combination 
probe. However, the characteristics of the probe in 
swirling flow were not fully determined. It was found in 
analyzing the results of the tests in the swirl annulus 
that the radial gradient in the flow angle gave rise to 
probe readings which were in error when the probe was ro- 
tated to balance. This was explained by the difference in 
the radial position of the sensors on the two faces of the 
wedge. It was therefore necessary to use the combination 
probe to sense flow direction in the compressor measure- 
ments and to set the wedge to the measured angle. 
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Table B1 

Probe Calibration Coefficients for =F0(X,0) 
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Table B2 



Error Analysis of Limiting 
Velocity From Combination Probe Calibration 
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SIDE VIEW 



FIG. B1 FREE-JET APPARATUS USED TO CALIBR,\TE PROBES 
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INCHES DOWNSTREAM, PROBE FROM RIGHT SIDE 
■ • M U M left " 
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NDN-D1MEN5IQNRL VELOCITY (X) 

FIG. B2 VELOCITY PROFILES FOR VARIOUS MACH NUMBERS MEASURED IN A 4.2 INCH FREE-JET 
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FIG. B3 ERROR IN PITCH ANGLE DERIVED FROM THE COMBINATION PROBE CALIBRATION 



(0) X=.I4B3MH 
(«) X=.lffi07fi 
(+) X=.2227B7 




.'iDN-i>!MEN5inNRL YEinClTY (X) 
(UHCDRREaED FDR SDUMDflRY LAYER EFFECTS) 
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(D) X=.IHg03£ RT 9.5 DIR. 

(*) X=.IH5I99 • • • 

(!) X=.22IBS • • • 
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NUMBER UF DIRMETER5 FROM PUSITIDN WHERE ERROR IN X=0 TO 
THE HRLL BDUNDRRY IN 10 INCH PIPE 

FIG. B5 ERROR IN VELOCITY DUE TO BOUNDARY EFFECTS IN A 10 INCH PIPE FLOW 




o 

M 

e 

o 

w 

a: 

a: 

o 

a 



o 

w 

/^rr, 

XCt, 

wpq 



— << 
«^Q 
az 

-J tZ) 

>PQ 



30 

C3P*h 



Sen 
! s 



ui 

s 



I 



m 

M 

O 






w 

o 

o 



VO 

PQ 



O 

HH 

(j^ 



I 



97 




98 



FROM BOUNDARY EFFECT 
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H.Bal 




ERROR IN PHI (PITCH RNBLE) 
DEE. 
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* DIAMETERS FROM HALL BOUNDARY IN IB INCH PIPE FLOW 
FIG. B8 ERROR IN PITCH ANGLE DUE TO BOUNDARY EFFECT IN A 10 
INCH PIPE FLOW 
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F" I m . 10 F3tt:QlJL.T£a OF* SH I HU. MMNIUU-UQ XGu£3X OF* COMM . HHdfBB: 



APPENDIX Ct INLET FLOW FIELD DETERMINATION 



Cl. METHOD OF APPROACH 

An important consideration in a transonic compressor 
installation is the velocity profile at the face of the 
rotor. The profile depends on the inlet shape and the mass 
flow rate. The following analysis is based on the assump- 
tion that a close approximation to the true velocity may 
be calculated if the referred flow rate is known. 

Seven surveys were conducted by Anderson (Ref, 3) in 
which accurate calculation of the inlet velocity profiles 
was accomplished using a calibrated United Sensor Probe 
over a range of flow rates tabulated in Table Cl, 

Reference 11 defines a "total flow function," , given 

by 

T = — ^ = F(X) = X(1 - C(l) 

^ ^t'^t U^t 

The flow function provides the means of replacing density 
in the equation with the non-dimensional velocity (X) for 
internal compressible flows. For each survey conducted in 
Hef. 3,"^ was calculated at each point and plotted versus 
the nonnalized displacement (hub to tip). The results are 
shown in Fig. Cl. Each curve was then fitted with a suit- 
able polynomial of the form 

^ = aQ + a^H + + a^H^ + ... C(2) 

where H = normalized displacement (hub to tip) and a^ , ag,.. 
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are determined by a polynomial regression curve fitting 
routine. The value of the total flow function at the cen- 



terline was evaluated for each case. Values of 

^/^CL for each survey were then calculated and plotted 
against H. The resulting curves shown in Fig. C2 were then 
approximated by a family of polynomials having the form 



1 =% 



c^ = Aq + A^H + A2H .... 

X / 

which is then called the normalized flow function. 



C(3) 



The referred flow rate is defined as 



where 



and 



^ dJVe 









S' 

"^ref 



■to 



rei 






and O = mass flow rate. 

A plot of^ versus '/if” shown in Fig. 
following linear relationship: 



C(4) 

C(5) 



C(6) 



C3 yielded the 



Thus, Eq. C(3) and Eq. C(?) provide an analytical rep- 
resentation of the inlet flow profile in terms of the 
referred flow rate only. 
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C2. APPLICATION 

, * 

For each run of the transonic compressor, W is calcu- 
lated and^^T is then evaluated from Eq, C(7). ^ can then 

CL 

be evaluated from Eq. C(3). 

X is then evaluated as the subsonic solution of the 
flow function given by Eq, C(l) using Newton's method as 
follows i 

1) Guess X = Xq 

2) Calculate 



^2 = X (1 - 


C(8) 


_ -ir 1 2 / ^ )] 

dX ^2 Lx ■ (y-1) ^ i_x2J 


C(9) 


X = X + 

dg 

dX 


C(10) 



3) Now test for convergence of 

If II 2 -5 00001 then the calculation is complete. 

If this requirement is not satisfied then the new I 2 
Eq. C(8) is calculated using the new X in Eq. C(10). 

4) Mach number is calculated . using Eq. B(4). 
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Table Cl 



Flow Rates Used for Inlet Velocity Determination 



Referred 
Flow Rate 


Centerline Value 
of Flow Function 


8.95 


.09589 


9.38 


.101607 


10.89 


.11693 


11.18 


.l'2039 


11.42 


.122025 


13.83 ' 


.151239 


14.30 


.156789 
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NI1N-01MEN5IDNBL DISTANCE (HUH TD TIP) 

FIG. Cl TOTAL FLOW FUNCTION VS. NORMALIZED DISPLACEMENT 
UPSTREAM OF THE ROTOR 



sri 
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NDHMflLIZED DISPLflCDIENT (HUH TD TIP) 

FIG. C2 NORMALIZED FLOW FUNCTION VS. NORMALIZED DISPLACEIVIENT 
UPSTREAM OF THE ROTOR 
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REFERRED FLDH RflTE 

FIG. C3 FLOW FUNCTION AT MID SPAN VS. REFERRED FLOW RATE 



APPENDIX Di ROTOR LOSSES 



Dl. NASA LOSS COEFFICIENT 

Losses incurred in flows in turbomachinery may be de- 
scribed with a varity of methods depending upon the measure- 
ments available and the specific application. Ref. 1 
describes a loss coefficient defined by 



tD’ = 







f 




D(l) 



where 

= mass-averaged defect in 
relative total pressure 
divided by the pressure 
equivalent of the inlet 
velocity head 

f 

P- = relative total pressure 
upstream 

Pp' = relative total pressure 
downstream 



P^ = static pressure upstream 

I 

(Pg)^ j = ideal relative total 
pressure downstream 




From the measurements obtained and a knowledge of the 



inlet velocity profile the parameters P^, P^, ^t2’ ^1* 

X T and T “ ' 

2 * ti' t2 are known. A representation of bJ in terms 

of measured stagnation conditions and velocities (X or M) 



is desired. 
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D1 , 1 Ups'tream of ~the Rotor . The following relations 
may be written! 






= P 



tl 




V^-1) 




y 2 

U2 J 



D(2) 



T 





D(3) 




Tti 




2 

+ X..,, 

W1 





D(4) 



where X is a velocity non-dimensionalized as in Sq, B(l), 
subscript denotes the relative velocity, and subscript W 
denotes the wheel speed. The other notation is given in 
the T-S diagram for the flow through the rotor shown in the 
above sketch. 

Using the above equations, for an "impulse" machine 
for which there is a meridional entrance flow, the following 
relations result! 



T 



E 



T 



tl 



1 + X 



2 

U2 



t 




Pti ( 1 + X 



2 

Ul 






D(5) 



D(6) 



and 

h = Pti ( 1 - x2 ) D(7) 

D1.2 Downstream of the Rotor . Sufficient measurements 

I 

exist at the machine station 2 to derive ^2 using 

^2*^2’ ^2 ^2 measured downstream velocity, 

non-dimensionalized using the downstream limiting velocity 
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(V^^) iri Eq. B(l). Alternatively, since Tg= constant for 
the flow through the rotor (see diagram), then one of the 
separate measurements is not needed. Two methods were 
therefore used. Method 1 used all measurements at station 
2, Method 2 did not use the temperature rise at station 2, 
and thus served as a check on the accuracy of the temper- 
ature measurement, 

(i) Method I (all measurements at station 2) 






D(8) 



where 






is given by 






F ? '^‘^1 

^ = 1 - X„ (sr^) - 2X^ X 



2 '"U3V T^2 



Tt2 ^2 '^t2 

Using Sq. D(6), D(7), D(8), and D(9) in Eq, D(l), there 

results 



D(9) 



Uj' 



[i.xy - 


r 


2X2Xu2SinC<2‘ 


^ — ) 


(Tti/Tta) ■ Tta/Tti J 








[MV-I ) . [1 . xf ] 


Y/(^-i) 



■(ii) Method II (Temperature measurement at station 2 

not needed) 

Temperature measurement downstream may be found using 
the follov/ing relation. 



•t2 



T^ 






Define 



. j. X 1C *t 2 ! — — — 

^ = XgXu 2 sinO^ sin 2 )^ + 1 



D(ll) 

D(12) 
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Using Eq. D(6), D(7)i D(ll), and D(12) in Eq, D(l), the 
loss coefficient is given by 



[i«y V(v-i). 

' [i + -[i - xf] 

D2.. DIFFUSION FACTOR 



D(13) 



Ref, 1 reports results of correlating compressor losses 
in teras of the "diffusion factor". This approach results 
from the argument that the loss in stagnation pressure de- 
pends largely on the suction-surface adverse pressure gra- 
dient, of which the diffusion factor can be shown to be an 
approximate measiore. The diffusion factor of Ref, 1 is 
defined by 



D = 1 



= — 4. 

< 5 ' (R^ + R^)W^ 



D(14) 



In this form, changes in axial velocity and radius through 
the blade row are accounted for. 

It is concluded in Ref, 1 that the total-pressure-loss 
parameter — ) can be correlated in terms of the 

above diffusion factor. Each family of blade shapes has a 
somewhat different correlation ciirve, however, this pre- 
sents no real problem as the total number of shapes is 
limited. The correlation for several types of blading is 
shown in Fig, 203 of Ref, 1, 
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